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Abstract 

Monte Carlo techniques have been widely employed in statistical physics as well as 
in quantum theory in the Lagrangian formulation. However, in some areas of appli- 
cation to quantum theories computational progress has been slow. Here we present 
a recently developed approach: the Monte Carlo Hamiltonian method, designed to 
overcome the difficulties of the conventional approach. 



1 Introduction 



The Monte Carlo (MC) method with importance sampling is an excellent 
technique to study systems with many degrees of freedom. It has been widely 
applied to statistical physics and quantum theories. There are two formula- 
tions of quantum theory: the Hamiltonian formulation and the Lagrangian 
formulation. A comparison of the conventional approaches is given in Tab. 
1. A natural question is whether one can perform efficiently MC simulations 
combined with the Hamiltonian formulation. If so, the disadvantages of the 
Hamiltonian formulation might be overcome. Our motivation is to construct 
an effective Hamiltonian via MC simulation starting from the original action. 
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Table 1 

Comparison of Hamiltonian and Lagrangian formulations. 



Formulation 


Hamiltonian 


Lagrangian 


Approach 


Solving Schrodinger Eq. 

TTl IP \ ZT" 1 ^ 

-n l-C/n >= -C'nI-C'n > 


Computing Path Integral 

^ /n ^_ / d[x]0[x]eM-Six]/K) 
J[dx]exp{-S[x]/n) 


Algorithm 


variational, series expansion, 


Monte Carlo simulation 




coupled cluster expansion ... 


as in Stat. Phys. 


Advantage 


can obtain not only energies 


use computer to generate 




of ground &: excited states, 


the most important configs. 




but also wavefunctions 


obeying Boltzmann's law 


Disadvantage 


methods above are too 


difficult to obtain: wavefunction. 




tedious for realistic 


excited states, S matrix. 




applications 


finite density QCD 



2 Algorithm 



Let us review briefly the basic ideas [1]. According to Feynman's path integral 
approach to quantum mechanics, the (imaginary time) transition amplitude 
between an initial state at position Xi, and time ti, and final state at Xf, tf is 
related to the Hamiltonian H by 

< Xf,tf\xi,ti >=< Xf\e~^^'^^~^'^^^\xi > 

oo 

= ^ < Xf\E, > e-"""^!^ < E,\xi >, (1) 
i/=i 

where T = tf — ti. The starting point of our method, as described in more 
detail in [1] is to construct an effective Hamiltonian i?e// (finite NxN matrix) 
by 

< Xf\e-"^^f-^'^/^\xi >^< Xf\e-"-ff^/^\xi > 

= ^ < XflEt^f > e-^S''^/^ < Eff\xi > . (2) 
i/=i 

Hgff can be found by MC simulation using the following procedure: 

(a) Discretize the continuous time. 

(b) Generate configurations [x] obeying the Boltzmann distribution 
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(3) 



(c) Calculate the transition matrix elements 



Mfi=< Xf\e-^^'ff^/^'\xi> (4) 

between N discrete Xi points and N Xf points. Note that the matrix M is 
symmetric. 

(d) Diagonalize M by a unitary transformation 



M^U^'DU, (5) 

where D = diag{e^^^^^'^/^'' , ...^e^^^^'^/^y Steps (a) and (b) arc the same as 
the standard MC method. Step (c) is the essential ingredient of our method, 
from which we can construct i^e//, and obtain the eigenvalues E^,^^ and wave- 
function \Elf^ > through step (d). Once the spectrum and wave functions are 
available, all physical information can be retrieved. Since the theory described 
by H is now approximated by a theory described by a finite matrix H^ff, the 
physics of H and -f^e// might be quite different at high energy. Therefore we 
expect that we can only reproduce the low energy physics of the system. This 
is good enough for our purpose. 



3 Testing the Method 



We have tested a number of quantum mechanical models in both 1+1 dimen- 
sions [1-3] and 2+1 dimensions [4,5]. Here we will show some new results for 
a model in 1+1 dimensions with the following potential 

V(x) = x^/2 + x^A. (6) 



The Euclidean action is given by S" = dt[mx^ /2 + This model is not 

exactly solvable. We have compared the results from our method with those 
from a standard algorithm (Runge-Kutta algorithm plus node theorem). The 
first three eigenvalues are given in Tab. 2. The parameters are m = 1, T = 1, 
h — 1, Ax — 1, N — 21. The results for the first two wave functions are shown 
in Figs. 1 and 2. Besides the eigenvalues El^^ and wave functions \El^f >, 
we have computed thermodynamical quantities such as the partition function 
Z, average energy E and specific heat C . Since we have approximated H by 
i^e//; we can express those thermodynamical observables via the eigenvalues 
of the effective Hamiltonian 
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Table 2 

Spectrum of the model with V{x) = a;^/2 + a;^/4. -E^ *' corresponds to Runge-Kutta 
plus node theorem and AE^-^- is the statistical error. 
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Tpn.t. 


T^M.C. 
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2 


3.6985 


3.6060 
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Fig. 1. Ground state wave function of the model with V{x) = x^/2 -|- x^/A. 



N 



v=l 



ii/=i 



(7) 



where (3 = T/h. The results for the average energy and specific heat are 
shown in Figs. 3 and 4, where the sohd hnes corresponds to the standard 
algorithm. Within statistical errors, our results are consistent with those from 
the standard algorithm. 



4 




Fig. 2. First excited state of the model with V{x) = jl + 




/3 



Fig. 3. Average Energy of the model with V{x) = + x'^/A:. 
4 Summciry and Outlook 



We have presented the basic ideas and some new results from our recently pro- 
posed Monte Carlo effective Hamiltonian method. We found that the method 
works well also for the V x'^ potential, allowing to compute excited states 
and thermodynamical obscrvablcs. In the previous simulations, the initial and 
final positions Xi and Xf were chosen to be uniformly distributed. As can be 
seen, for a given dimension of the Hilbcrt space N , the results break down if j3 
is too small. Increasing N requires larger CPU time. For quantum mechanics 
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0123456789 10 
Fig. 4. Specific heat of the model with V{x) = + 



in 3+1 dimensions or quantum field theories, a stochastic basis will be neces- 
sary to select the most important contributions to the transition matrix Mjj. 
Such a work is in progress [3]. 
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